Attosecond control of electron dynamics in carbon monoxide 
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Laser pulses with stable electric field waveforms establish the opportunity to achieve coherent 
control on attosecond timescales. We present experimental and theoretical results on the steering of 
electronic motion in a multi-electron system. A very high degree of light-waveform control over the 
directional emission of and fragments from the dissociative ionization of CO was observed. 
Ab initio based model calculations reveal contributions to the control related to the ionization and 
laser-induced population transfer between excited electronic states of CO^ during dissociation. 

PACS numbers: 



Coherent control of chemical reactions and photobio- 
logical processes has been achieved by manipulating the 
laser frequency, phase and polarization in closed loop 
experiments [1]. Control of the electric field waveform 
E{t) = Eo{t) cos{ujt + (/)), with envelope Eo{t), and fre- 
quency cc;, by the carrier envelope phase (CEP) (j) consti- 
tutes a new paradigm of coherent control. This control 
became accessible with CEP-stabilization and opened 
the door to steer electrons in atomic and molecular 
systems on attosecond timescales [2]. Waveform con- 
trolled few-cycle pulses have only recently been used to 
control electron localization in the dissociative ioniza- 
tion of the prototype molecules D2 (1, H] and HD 
The same processes were also theoretically studied (see 
e.g. ^51, |6|, i2, i8|, i^] and references cited therein). After ini- 
tial ionization, these systems only contain a single elec- 
tron. The important question arises, whether the steer- 
ing of electrons in more complex systems is feasible and 
— if yes — can we understand the role of the initial ion- 
ization/excitation process and the following strong-field 
coupling of the various potential energy surfaces in the 
observed control? We describe experiments, where con- 
trol of electron dynamics in carbon monoxide (CO) was 
achieved by the light- waveform. Phase-stabilized 4 fs, lin- 
early polarized laser pulses at 740 nm and at an intensity 
of 8 X 10^^ W cm~^ were applied to dissociatively ion- 
ize CO. The directional emission of ionic fragments was 
monitored via velocity-map imaging ( VMI) . We compare 
the experimental results to full quantum calculations al- 
lowing a mechanistic interpretation and understanding 
of the observed control. The dynamics of molecules in 
strong laser fields typically includes ionization and dis- 
sociation. While there is a wealth of work on the dis- 
sociation of small molecules in strong laser fields (see 
e.g. [10] and references cited therein), only a limited num- 
ber of studies have been performed on CO. Guo studied 
the multiphoton induced dissociative ionization of CO 
at 800 nm and an intensity of 3.8 x 10^^ W cm-^ 
The yield of C+ ions from the dissociation was found to 
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Figure 1: PES of C0+ and Frank Condon point (FC) for the 
ionization from neutral CO obtained by calculations described 
in the text. 



be more than an order of magnitude higher than that 
of 0+ fragments. Experiments reported by Alnaser et 
al. |l2| using 8fs pulses at 800 nm and an intensity of 
6 X 10^^ W cm~^ showed an angular distribution of the 
ionic fragments from the Coulomb explosion of CO with a 
maximum along the laser polarization axis, for which the 
shape of the highest occupied molecular orbital (HOMO) 
of the CO molecule was held to be responsible. Here, 
we report on the CEP-dependent directional emission of 
ionic fragments from the dissociative ionization of CO. 
To demonstrate the complexity of the molecular system, 
we show in Fig.[T]the calculated potential energy surfaces 
(PES) of the lowest ^E+- and ^H-states relevant for the 
photodissociation of C0+ at our laser intensities. In the 
initial ionization, the low-lying bound electronic states X 
and A of C0+ can be populated. Higher electronic states 
of C0+ can be reached by recollision or multi-photon ex- 
citation starting the dissociation process. At our inten- 
sity the recollision energy of the electron from the ini- 
tial ionization is up to 13 eV. The setup used in these 
experiments was described earlier Briefiy, phase- 

stabilized few-cycle pulses were focused into the center 
of the ion optics of a VMI spectrometer using a spherical 
mirror (f = 125 mm). Ions and electrons that were gener- 
ated at the crossing point of the laser (linearly polarized 



2 



along the y-axis and propagating along the x-axis) and 
an effusive atomic/molecular beam were projected (along 
the z-axis) onto a multi-channel plate (MCP) /phosphor 
screen assembly and recorded with a cooled CCD cam- 
era. Inversion of the recorded projections using an it- 
erative procedure allowed reconstruction of the original 
3D ion momentum distributions. Fig. [2](a) shows a cut 
through the 3D momentum distribution in the x^-plane 
at Pz=0 for C+ ions and for a CEP:^ tt pulse. An up- 
down (positive py versus negative py values) asymmetry 
in the C+ fragment emission is visible and the contribu- 
tions may be best identified as three rings, where the first 
ring is the broadest (from p = to l-10~^^A^s) and most 
intense and exhibits additional sharp lines. The other 
rings appear between momenta of 1.1 to 1.3 and 1.4 to 
1.6-10~^^iV5. The dominant features show angular distri- 
butions with minima around 0°, 90°, 180° and 270° with 
varying modulation depths. This angular distribution of 



the C+ fragments differs from the findings in ref. [12| 
and indicates contributions from the HOMO (5cr) and 
the HOMO-1 (iTT^/y). The solid black line in Fig. E^a) 
displays the calculated volume-averaged angle-dependent 
ionization rate as described below in the theory section. 
Fig. [2Kb) shows the kinetic energy spectrum derived from 
Fig. EJa) by angular integration (blue line) together with 
the calculated spectrum obtained within the theoretical 
approach described below (red line). Comparison with 
experiments using circularly polarized light at twice the 
intensity as for linear polarization, which suppresses the 
recollision excitation, reveals a distinct decrease in in- 
tensity of all three rings. Thus, we conclude that rec- 
ollision excitation is responsible for the production of 
the observed ionic fragments at our experimental con- 
ditions. In the fragmentation of C0+, the C++0 chan- 
nel is energetically favored over the C+0+ channel (see 
Fig. [1]). We find the C+ yield to be approx. 20 times 
larger than the 0+ yield. Throughout the measured 
CEP range a pronounced phase dependence in the direc- 
tional ion emission is found for C+ as well as for 0+ ions. 
The directional emission is represented by the asymmetry 

ment kinetic energy W and the laser phase (j). Pup{W^ (j)) 
and Pdown(^, (/>) are the angle-integrated ion yields in 
the up and down directions. Fig. [2Kc) displays the ob- 
served asymmetry A{W^ (j)) for the dissociative ionization 
of CO into C+ and O as a function of the CEP and the 
kinetic energy W of the C+ ion fragments. Fig. [2](d) 
shows the asymmetry parameter integrated over selected 
energy ranges. The observed asymmetry in the direc- 
tional emission of C+ ions is very pronounced and al- 
most equally strong throughout the kinetic energy spec- 
trum. A similar asymmetry map (not shown here) was 
recorded for 0+ ions, showing the same features. Due to 
the significantly weaker 0+ signal, however, the asym- 
metry map exhibits a lower signal-to-noise ratio, making 
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Figure 2: (a) Inverted two-dimensional momentum image 
(the laser polarization is vertical) and calculated ionization 
rate of CO in arb. un. (black line) for CEP~ tt; (b) measured 
(blue) and theoretical (red) C^ kinetic energy spectrum for 
CO dissociative ionization, (c) Asymmetry of C^ ion emission 
along the laser polarization (integrated over 120°) vs. kinetic 
energy and phase, (d) Asymmetry integrated over indicated 
energy ranges vs. CEP. The CEP was calibrated by reference 
measurements in Xe and comparison to calculations based on 
the quantitative rescattering theory according to a recently 
published method [l^l. From the evaluation of this data we 
obtain an error in the absolute CEP of di0.047r. 



it also difficult to determine the phase-shift between the 
points, where a maximum asymmetry is found for C+ and 
0+ ions. In the experiment, the asymmetry can arise 
from contributions of all three steps: ionization, recol- 
lisional excitation and laser-induced population transfer 
between excited electronic states of C0+. The initial 
population of excited states by recollision is likely de- 
pendent on the CEP, but its calculation is currently out 
of scope for larger molecules @. We can, however, focus 
on the remaining two steps. The ionization probability 
of a molecule in a laser field is determined by the electron 
fiux induced by this external electric field [l5|. To cal- 
culate the angular dependent ionization probability for a 
given static electric field, we performed quantum chemi- 
cal calculations under various orientation angles with re- 
spect to the polarization of the applied external field. 
In the spirit of [15] we record the electron fiux through 
a surface placed at the outer turningpoints of a given 
orbital. When ionization from the HOMO only is con- 
sidered as in ref. (l2| . our calculated angular distribution 
matches the reported results therein. For our experi- 
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mental conditions we allow ionization from the HOMO 
(leading to the X state of C0+) and the HOMO- 
1 (leading to the A ^E-state) reproducing the observed 
angular distribution (Fig.[2](a)). This is achieved by a ba- 
sis transformation forming the orbitals HOMO+HOMO- 
1 and HOMO-HOMO-1. By integration of the calculated 
ionization probabilities over the appropriate orientation 
angles, the asymmetry in the ionization step caused by a 
phase-stabilized electric laser field can be extracted. As- 
suming ionization only at the extrema of the field, we 
calculated the asymmetries for various contributions of 
half-cycles as shown in Fig. [3] (solid line). The volume- 
averaged angular dependent ionization probability is ob- 
tained for a spatial gaussian intensity profile and is shown 
in Fig. [2](a) (black line) for the combination of the three 
most prominent half-cycles and a CEP=7r pulse. The 
related volume-averaged asymmetries are presented in 
Fig. [3] (dashed line). The remaining less intense half- 
cycles can be neglected as they do not cause recollision 
excitation. The good agreement between theory and ex- 
periment shows that the angular distribution arises from 
the ionization out of two orbitals. Its asymmetry origi- 
nates from phase stable electric fields. The corresponding 
calculated asymmetry amplitude of 0.14 is below the ex- 
perimentally observed total asymmetry amplitude of 0.2. 
If only ionization would be responsible for the observed 
CEP-dependence, the asymmetry would be expected to 
be seen for all observed C+ ions irrespective of their ki- 
netic energy or angular distribution. Experimentally, we 
do however find that the angular range over which the 
asymmetry is observed depends on the fragment kinetic 
energy, suggesting that the experimental result is not 
solely explained by the ionization mechanism. To calcu- 
late to what extent the asymmetry may arise in the disso- 
ciation process we used our recently introduced approach 
to describe coupled electron and nuclear dynamics 
We follow the dominant recollision excitation pathway, 
assuming ionization of CO only at the electric field maxi- 
mum of the laser pulse used in the experiment and recolli- 
sion of the electron after 1.7 fs [16]. Our quantum dynam- 
ical simulations start after the recollision. The initially 
formed nuclear wave packet is excited from the X and 
A state to higher electronic states. Thereby a coherent 
superposition of several electronic states is created and si- 
multaneously the electronic and nuclear wavepacket mo- 
tion is initiated. The PES for the CO electronic ground 
state and for the ionic states are calculated with MOL- 
PRO [17]. From the set of C0+ PES we chose three poten- 
tials as typical representatives for the induced dynamics 
which moreover allow all transitions among each other. 
We include the C state as the weakly bound state 
and the E ^11 state to resemble the repulsive dynamics. 
Both states correlate with the C+(^P)+0(^D) channel. 
As third PES we include the H state which is the 
first state leading to the C(^P)+0+(^S) channel, deliv- 
ering O^ fragments as observed in the experiment. The 
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Figure 3: Calculated asymmetries (+ symbols and solid line) 
for various combinations of half-cycles of the laser field (shown 
in the inset). The x symbols show volume averaged asymme- 
try contributions (dashed line). 
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Figure 4: (a) Electric field (CEP = (p=0). (b) Time- 
dependent populations of the C E and H states 
after recollision excitation (solid: CEP=0; dotted: CEP=7r). 
(c) Temporal evolution of the probability measuring a C^ 
fragment Pc+- 



initial wavepacket is composed as a 55:38:7 distribution 
of the states involved, ordered in increasing energy, as- 
suming a gaussian energy distribution for the recollid- 
ing electron with a cut off energy of 13 eV. Choosing 
these particular values for the initial populations also 
reflects the experimentally observed ratio between C+ 
and 0+ fragments. In the calculations the C0+ ions are 
taken to be aligned at 45° to the laser polarization as 
the ionization peaks along this direction and all transi- 
tions between S and 11 states are allowed. The molecular 
wavefunction ^moi = '^iXi{Rit)(pi{r^t] R) is set up as 
the sum over these three electronic states i with Xi the 
nuclear wavefunctions, cpi the fleld-free electronic wave- 
functions, the nuclear and electronic coordinates R and 
r and the time t. Figure [4] shows the temporal evolution 
of the laser fleld (a) and of the population in the three 
selected electronic states (b). The calculated kinetic en- 
ergy spectrum (Fig. Efb)) derived from the nuclear dy- 
namics is in reasonable qualitative agreement with the 
experimental data supporting the selection made for the 
representative states. Note that contributions from dou- 
bly charged states cannot be completely ruled out in the 
experimental spectrum. On the basis of these results, we 
can explain the origin of the observed energy distribu- 
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tion. The low kinetic energy spectrum arises from the 
dynamics on the weakly bound state while the spectrum 
between 2.0-2.8 eV originates from the purely repulsive 
state. The high energy spectrum reflects the dynamics of 
the second repulsive state correlating with the 0+ chan- 
nel. The sharp peaks in the low kinetic energy spectrum, 
also present in the experimental data (Fig. [2](a,b)), arise 
from interference of bound and dissociative vibrational 
states in the nuclear wave packet on the C state. 
The electronic density p{ri , t; R) is expressed as a func- 
tion of the electronic coordinate ri and time. The total 
molecular wavefunction is integrated over the nuclear and 
over the — 1 electronic coordinates (with N the total 
number of electrons): 

3 

= ^a,{t)^\^,{n,to;R)\^ 

(1) 

3 ^ ^ 

-^^^'^Re[{x^{R^t)\xJ{R^t)), 

i—1 j<i 

^,(ri,to;i^)^,(ri,to;i?)e-^^^^-^(^(^))^W-At)| 



with ai{t) = y {xi{t)\Xi{t)) R and af(t) the popu- 
lation of the electronic states, the interference term 
{Xi{t)\xj{t)) R and the energy difference AEij{t) = 
Ej{t) — Ei{t) between the electronic states i and j. The 
electronic wavefunctions (pi{r^to;R) are represented as 
Slater determinants. The time evolution of the elec- 
tronic wavepacket is calculated by propagation in the 
eigenstate basis. The coupling of the fast electron to 
the slower nuclear dynamics enters through the time- 
dependent population and through the interference term, 
which specifles the degree of electronic coherence in- 
duced in the molecular system. The probability Pc+{t) 
of measuring a C+ fragment is given by Pc+{t) = 
/;-^;dx/;--;d^/^^^^^ where x,y and z 

refer to the molecular frame with z along the molecular 
axis and the 0-atom oriented along negative values 
In the flrst 6 fs the electron dynamics, reflected in Pc+ (t) 
(e.g. red curve in flg. ID^c)), results from a competition 
between the influence of the light pulse, the dynamics of 
the linear combination and the interference term of the 
nuclear wavefunctions. As soon as the population trans- 
fer between the electronic states stops (approximately 
after 8 fs, see flg. SKb)) Pc+{t) oscillates with decreas- 
ing amplitude converging after 12 fs to its flnal value. 
This decay arises from the reduced overlap of the super- 
imposed molecular orbit als, which become soon located 
on the two different nuclei during the dissociation pro- 
cess. Consequently, the last term of eq. [T] vanishes and 
hence the dynamics of the electronic linear combination. 
The probability Pc+{t) is flnally given by the prepared 



population distribution of the coupled reaction channels 
leading either to C+ or 0+. This ratio is steered very 
precisely by the CEP modulating the relative intensities 
between half-cycles and is only signiflcant in ultrashort 
laser pulses. When multiple half-cycles are taken into 
account the mechanism of the observed control stays the 
same. A shift of the CEP by tt while keeping the molec- 
ular orientation leads to a different result (green curve in 
flg. SJc)). Changing the orientation of the molecule by 
180° is equal to shifting the CEP by tt as the transition 
dipole moment changes the sign. Thus the CEP depen- 
dent asymmetry in the dissociation step can be calculated 
by Pc+{t) for two CEP values shifted by tt. As shown 
in flg. SKc), the probability of measuring a C+ fragment 
upon the break up of the molecule can be changed by 5% 
through the CEP. 

We presented experimental and theoretical results on 
the steering of electrons in a multi-electron system. As 
possible mechanisms for the observed CEP-control of the 
directional fragment emission in the dissociative ioniza- 
tion of CO we have discussed contributions from the ion- 
ization step as well as from the laser induced dynamics 
during the dissociation. The current experimental data 
does not allow to clearly distinguish the individual con- 
tributions. Further studies are underway that aim into 
this direction. 
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